Group-theoretical construction of extended baryon operators in lattice QCD 
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The design and implementation of large sets of spatially-extended, gauge-invariant operators for 
use in determining the spectrum of baryons in lattice QCD computations are described. Group- 
theoretical projections onto the irreducible representations of the symmetry group of a cubic spatial 
lattice are used in all isospin channels. The operators are constructed to maximize overlaps with 
the low- lying states of interest, while minimizing the number of sources needed in computing the 
required quark propagators. Issues related to the identification of the spin quantum numbers of the 
states in the continuum limit are addressed. 
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I. INTRODUCTION 

Spectroscopy is a powerful tool for uncovering the im- 
portant degrees of freedom of a physical system and the 
interaction forces between them. The spectrum of quan- 
tum chromodynamics (QCD) is indeed very rich: con- 
ventional baryons (nucleons, A, A, S, f2, . . .) and mesons 
(•7T, K, p,4>, . . .) have been known for nearly half a cen- 
tury, but other higher-lying exotic states, such as glue- 
balls, four-quark states, and so-called hybrid mesons and 
hybrid baryons bound by an excited gluon field, have 
proved more elusive, mainly because our theoretical un- 
derstanding of such states is insufficient, making their 
identification problematical. 

Recent discoveries of several new hadronic resonances 
have generated much excitement in the field of hadron 
spectroscopy. The E852 collaboration^ reported a sig- 
nal for a 1 ^ hybrid meson decaying into pir with a mass 
near 1.6 GeV, though the significance of this result has 
been questioned 0. Another exotic 1 ^ candidate at 1.4 
GeV has been tentatively identified in the ijir channel 
by E852 3], VES4], KEK 5], and Crystal Barrel 6]. The 
first observation of a doubly-charmed baryon has been 
reported 7], and evidence for the possible existence of 
a strangeness 5 = 1 pentaquark state has emerged at 
SPring-8, JLab, and elsewhere |sl [Tol FTT| ■ Interest in 
excited baryon resonances has also been raised by experi- 
ments dedicated to mapping out the N* spectrum in Hall 
B at the Thomas Jefferson National Accelerator Facility 
(JLab), and the search for hybrid mesons and glueballs 
is intensifying due to the Hall D initiative at JLab and 
experiments at CLEO-c. 

Unfortunately, our understanding of many conven- 
tional and excited hadron resonances is vestigial and 
comes only from QCD-inspired phenomenological mod- 
els, such as the bag model, the nonrelativistic quark 
model, and quark-diquark models, or from approaches, 
such as QCD sum rules and methods based on Schwinger- 
Dyson equations, which use approximations whose justi- 



fications are unclear. There are a growing number of 
resonances which cannot be easily accommodated within 
quark models. States bound by an excited gluon field, 
such as hybrid mesons and baryons, are still poorly un- 
derstood. The natures of the Roper resonance and the 
A(1405) remain controversial. Experiment shows that 
the first excited positive-parity spin- 1/2 baryon lies be- 
low the lowest-lying negative-parity spin-1/2 resonance, a 
fact which is difficult to reconcile in quark models. Given 
the current surge in experimental activity, the need for an 
understanding of such states from QCD itself has never 
been greater. 

This need has motivated us to undertake a compre- 
hensive ab initio study of the hadron spectrum in QCD. 
Presently, Monte Carlo estimates of QCD path integrals 
defined on a space-time lattice offer the best way to make 
progress in this regard, so this is the calculational ap- 
proach we have adopted. The Monte Carlo approach has 
been used to investigate hadrons throughout the past 
two decades, but the number of states studied to date 
has been somewhat limited. Also, prior works have not 
strived to identify the continuum spin J of the states 
studied, simply assuming that the lowest allowed J would 
be the lowest-lying state. The novel features of our ap- 
proach are its comprehensiveness and its use of tech- 
niques to identify spin. Given the vast amount of ex- 
perimental data being generated at accelerator facilities, 
such as Jefferson Lab, there is an urgency to investigate 
the spectrum (masses, widths, transition rates, form fac- 
tors, and so on) as completely as possible. The number 
of hadron eigenstates which can be reliably extracted in 
Monte Carlo computations is not currently known, so our 
undertaking will be partly an exploration of the limits of 
possibility. Another aim is to discover whether a very 
large number of interpolating operators are needed to 
extract the low-lying spectrum, or whether a handful of 
carefully chosen ones is sufficient, and this work outlines 
a systematic means of finding such operators. 

Our first goal is to calculate the masses of as many 
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low-lying hadron resonances as possible in QCD using 
Monte Carlo techniques. The masses and widths of reso- 
nances (unstable hadrons) cannot be calculated directly 
in finite- volume Monte Carlo computations, but must be 
deduced from the discrete spectrum of finite- volume sta- 
tionary states for a range of box sizes |l2lll3lll4lll5| . The 
rigorous application of such techniques to obtain the res- 
onance parameters to high accuracy would require vast 
computational resources, but our goal here is merely to 
obtain a first exploratory scan of the spectrum of QCD, 
not to pin down each mass to very high precision. Hence, 
simply obtaining the finite-volume spectrum for a few 
judiciously-chosen volumes should suffice for ferreting out 
the hadron resonances from the less interesting scattering 
states and may even give qualitative information about 
preferred decay modes. 

To compute the finite- volume stationary state energies, 
the temporal correlations Cy(i) = (0\T O l (t)O j (0)\0), 
where T denotes time-ordering, of a set of operators 
Oj (0) which create the states of interest at an initial time 
t = with a corresponding set of operators Oi (t) which 
annihilate the states of interest at a later time t must 
be determined. The correlation functions Cy(i) can be 
expressed in terms of path integrals over the quark and 
gluon fields, and when formulated on a Euclidean space- 
time lattice, such path integrals can be estimated using 
the Monte Carlo method with Markov-chain importance 
sampling. Incorporating the quark-field effects into the 
Monte Carlo updating for realistically light quark masses 
remains a challenge, but there is steady progress with im- 
proving algorithms and increasing computational power. 

The procedure for extracting the lowest stationary- 
state energies Eq, E%, E%, . . . from the hermitian matrix 
of correlation functions Cij(t) is well known [TfillT^ I. Let 
A n (£, to) denote the eigenvalues of the hermitian matrix 
C{t )-^ 2 C(t) C{to)~ 1/2 , where t is some fixed reference 
time (typically small) and the eigenvalues, also known as 
the principal correlation functions, are ordered such that 
Ao > Ai > • • • as t becomes large. Then one can show 
that 

limA„(Mo) = e - £ " (t -* 0) (l + O( e - A " ( *- t0) )), (1) 
A„ = min\E k ~E n \. (2) 

Determinations of the principal correlators X n (t,t ) for 
sufficiently large temporal separations t yield the desired 
energies E n . 

The above equations illustrate the difficulties which 
must be faced in order to extract the stationary state 
energies Eq, E\, . . . from the temporal correlations of the 
hadronic operators. In each principal correlator, there 
are contaminating contributions from all other states 
which can be created and annihilated by the operators 
used. In order to reliably extract the single decaying ex- 
ponential of interest, the obscuring contributions from all 
of these other states must somehow be suppressed. 

There are two crucial ingredients in reducing the un- 
wanted contributions in \ n {t,to). The first is to use 



sufficiently large values of t. However, there are often 
practical considerations which limit how large t can be, 
and the statistical uncertainties in the Monte Carlo es- 
timates of the correlation functions generally increase 
with t. The second, and more important, consideration 
in suppressing the contamination in \ n (t, to) is to use 
cleverly- devised operators which couple minimally to the 
unwanted states. 

Successfully extracting the spectrum of QCD in our 
Monte Carlo computations will hinge crucially on using 
carefully designed hadronic operators. Excited meson 
and baryon resonances are expected to be large objects, 
so the use of spatially-extended operators is important. 
Since our calculations will be carried out on hypercubic 
space-time lattices, the energies of states in all irreducible 
representations (irreps) of the Oh cubic point group must 
be determined in order to identify the continuum-limit 
spin J of each physical state. Because determining the 
mass of a particular resonance requires determining the 
energies of all lower-lying stationary states, including 
scattering states, the set of operators we use must include 
not only meson and baryon operators, but also multi- 
hadron operators. Another very important fact to keep 
in mind is the computational cost of evaluating quark 
propagators, especially for light quark masses. Hence, 
our operators must be devised with an eye towards min- 
imizing the number of sources needed to calculate the 
required quark propagators. 

Designing the hadronic operators is an important first 
step in our comprehensive study of the mass spectrum 
in lattice QCD. From the above considerations, the guid- 
ing principles in devising our operators are maximizing 
overlaps with the states of interest while minimizing the 
number of quark-propagator sources. Although opera- 
tors for baryons, mesons, and their scattering states will 
be needed, we restrict our attention to the construction 
of the three-quark baryon operators in this first paper. 
Due to the complexity of these calculations and the im- 
portance of providing checks on our final results, we have 
been pursuing two different approaches to constructing 
the baryon interpolating field operators. Only one ap- 
proach is described here; an alternative approach based 
on Clebsch-Gordan techniques will appear elsewhere [l8|. 
Meson and multi-hadron operators will be detailed in 
subsequent works. Furthermore, this paper deals only 
with issues related to the construction and utilization of 
these operators. Results from Monte Carlo calculations 
using these operators will be presented in later publica- 
tions. 

This paper is organized as follows. An overview of our 
approach to constructing the hadronic operators is first 
outlined in Sec. [H] Our operators are assemblages of ba- 
sic building blocks, which are described in Sec. lIIII along 
with the conventions we use for the Dirac 7-matrices 
and in Wick rotating into Euclidean space-time. We use 
quark fields which are Dirac spinors, so our hadron op- 
erators apply to Monte Carlo calculations involving Wil- 
son, domain-wall, and overlap fcrmion actions, but not 
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to computations involving staggered fermions. The basic 
building blocks are then combined into gauge-invariant 
three-quark operators (referred to as elemental opera- 
tors) having appropriate flavor structure in Sec. If VI Pro- 
jections onto the rotation-reflection symmetry sectors 
produce the final operators in Sec. Issues related to 
the use of these projections in constructing the baryon 
propagators are discussed in Sec. IVII Concluding re- 
marks and plans for future work are outlined in Scc. lVIII 



II. OVERVIEW OF OPERATOR 
CONSTRUCTION 

Devising relativistic hadronic operators in continu- 
ous space-time usually involves combining Dirac spinors 
to form Lorentz scalars, pseudoscalars, vectors, axial- 
vectors, and so on, using the Dirac 7-matrices and the 
charge conjugation matrix C satisfying C 7^ C _1 = — 7J. 
It has been common practice in lattice QCD simulations 
to use hadron operators built in a similar fashion through 
simply discretizing their continuum analogs. However, 
such an approach becomes very cumbersome when con- 
structing higher spin states or complicated extended op- 
erators. Also, the above operators generally couple to 
states belonging to different J p (spin-parity) symmetry 
sectors. Since the hypercubic lattice breaks Lorentz co- 
variance, there is really no reason to construct operators 
according to Lorentz symmetries. Finally, we wish to ex- 
tract a large portion of the low-lying spectrum, which 
means that large sets of operators will be needed to com- 
pute complete correlation matrices. Hence, the usual ap- 
proach which mimics that used in continuous space-time 
is not feasible for our purposes. 

Instead, we advocate an approach which more directly 
combines the physical characteristics of baryons with the 
symmetries of the lattice regularization of QCD. Recall 
that baryon states are characterized by their total mo- 
mentum p, their total (half-integral) spin J, a projection 
A of this spin onto some axis (the z-axis or the momen- 
tum, say), their parity P = ±1, and their quark flavor 
content. The masses of the light u and d quarks are very 
nearly equal, and therefore, we work in the approxima- 
tion that m u = md- In this approximation, the theory 
has an exact isotopic spin symmetry, and states carry two 
more labels, total isospin / and its projection I3 onto a 
given axis. The other flavor quantum numbers which 
label the states are strangeness S, charm C, and bot- 
tomness B (we do not consider the top quark). In our 
calculations, isospin remains an exact symmetry since we 
neglect electromagnetic interactions. 

Since our simulations will be performed using a hyper- 
cubic space-time lattice, our operators should be classi- 
fied according to the symmetries of the lattice, rather 
than the full rotational symmetries of continuous space- 
time. Of course, we expect to recover the symmetries 
of the space-time continuum as the lattice spacing is 
made small. Since we are interested only in determin- 



ing the masses of the baryon states, we restrict our at- 
tention to representations corresponding to zero total 
three-momentum p = 0. Hence, our operators must 
be invariant under all allowed spatial translations and 
we require that they transform under spatial rotation- 
reflection symmetry operations according to the irre- 
ducible representations of the octahedral point group Oh- 
These irreducible representations are the lattice analogs 
of the continuum J p labels, and the row of the represen- 
tation is the analog of the spin projection A. Thus, our 
(annihilation) operators can be written 

Bt XF (t) = "£Bt XF (x,t), (3) 

X 

where A indicates the irreducible representation of Oh, A 
is the row of the A representation, F denotes all of the 
quantum numbers associated with the flavor content of 
the operator, and i labels the different operators in the 
AXF symmetry sector. Under a symmetry operation R, 
these operators transform according to 

U R Bt XF {t) 4 = E B i" F ® T< £( R )*> ( 4 ) 

A" 

u R Bt XF (t)ui = E^*) r$(R), (5) 

I" 

where Ur denotes the quantum operator which effects 
the symmetry operation corresponding to group element 
R (not to be confused with the gauge link variables), 
and r^(i?) are the elements of the A representation 
matrix corresponding to group element R. For baryons, 
we are only interested in states corresponding to half- 
integral spin J in the continuum limit, so we can restrict 
our attention to the six spinorial representations of Oh- 
There are four two-dimensional irreducible representa- 
tions Gig,Gi u ,G2g, and G2 U (adopting a Mullikcn-likc 
naming convention), and two four-dimensional represen- 
tations H g and H u . These representations will be dis- 
cussed in greater detail later. 

Our general approach to constructing the B^ XF (t) op- 
erators is to (1) first identify appropriate basic "building 
blocks" to use in constructing all baryon operators, (2) 
devise simple elemental operators containing the appro- 
priate flavor and color structure, then (3) apply appropri- 
ate group-theoretical projection operators to find linear 
combinations of the elemental operators with the desired 
transformation properties under the symmetry group of a 
spatial cubic lattice. Let Bf (t) = J2 X Bf (x, t) denote a 
gauge-invariant elemental operator with the appropriate 
quark flavor content and which is invariant under allowed 
spatial translations, then an operator which transforms 
according to the row A of the A irreducible representation 
is obtained using 

Bt XF (t) = — £ rf A >(i?) U R Bf{t) ul (6) 

where O® is the double group of Oh, R denotes an el- 
ement of Of, g D is the number of elements in Of , 
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and d\ is the dimension of the A irreducible represen- 
tation. Projections onto the double-valued irreps of Oh 
require using the double group Of L in Eq. ijfjj. Given 
Mb elemental Bf operators, many of the projections in 
Eq. © vanish or lead to linearly-dependent operators, 
so one must then choose suitable linear combinations of 
the projected operators to obtain a final set of indepen- 
dent baryon operators. Thus, in each symmetry channel, 
one ends up with a set of r operators given in terms of a 
linear superposition of the Mb elemental operators: 

M B 

Bt XF (t)=J2& F Bf(t), i = l...r. (7) 

Note that the expansion coefficients in the B t opera- 
tors are the complex conjugates of those in B^ XF : 

Mb 

Bf XF (t)=J2 ^ F *Bj(t), i=l...r. (8) 

III. THE BASIC BUILDING BLOCKS 

The oscillating path integral weight e lSM in quantum 
field theory, where Sm is the action defined in Minkowski 
space-time and using natural units h = c = 1, is unsuit- 
able for applying the Monte Carlo method to evaluate 
the correlation functions of the theory via Feynman path 
integrals. A rotation to imaginary time t — > —it, where r 
is real, leads to path integrals with weight e~ s , where S 
is the action defined in Euclidean space-time. If S is real, 
the path integral weight is real and positive and, hence, 
can be interpreted as a probability, allowing the applica- 
tion of Monte Carlo methods with suitable importance 
sampling. The Euclidean action S is defined such that S 
is invariant under all symmetries of Euclidean space-time 
and all Green's functions of the theory are identical to 
the Green's functions of the Minkowski theory, analyti- 
cally continued to imaginary time t — > —it. Although 
our simulations employ a path integral quantization of 
the field theory, a canonical quantization viewpoint can 
be adopted when discussing the quantum operators. 

Our conventions for the continuation from Minkowski 
space-time with metric g^ v — diag(l, — 1, — 1, — 1) into 
Euclidean space-time (imaginary time) are as follows. 
We define the following Euclidean space-time coordinates 
and derivatives (a subscript or superscript M indicates a 
Minkowski space-time quantity): 

X — X ^ — - IX ■ X^ — Xj — - X ji^j — X ' j . 

d i = d 4 = -id$", & = dj = -& M = df , (10) 

for spatial directions j — 1,2,3. The metric in Euclidean 
space-time is 8^ so there is no distinction between co- 
variant and contravariant indices. Our Monte Carlo cal- 
culations will be carried out using a hypercubic space- 
time lattice, and we require that the lattice spacings in 



the three spatial directions are the same, denoted by a s ; 
the temporal lattice spacing at may differ from a s , al- 
lowing us to exploit the known benefits of anisotropic 
lattices^). Throughout this paper, we set a s = 1 to 
simplify the notation. Four-vectors of unit length point- 
ing along the spatial axes of the lattice with a vanishing 
temporal component will be denoted by j, k, and so on, 
for j,k = ±1,±2,±3. 

As usual in lattice gauge theory, the gluon field is in- 
troduced using the parallel transporter J7 M (x) given by 
the path-ordered exponential of the gauge field along 
each link connecting neighboring sites of the hypercu- 
bic lattice. We also introduce the Dirac spinor field 
ipaai x ) which annihilates a quark and creates an anti- 
quark, where A refers to the quark flavor, a refers to 

A 

color, and a is the Dirac spin index, and the field ip aa (x) 
which annihilates an antiquark and creates a quark. Un- 
like in Minkowski space-time, ip and ip must be treated 
as independent fields, so we emphasize that tp ^ 
This is required in order to simultaneously satisfy Eu- 
clidean covariance of the fields, the canonical anticom- 
mutation relations, and the equality of the Euclidean 
two-point function with the relativistic Feynman prop- 
agator continued to imaginary times po). l2l| . Our Eu- 
clidean space-time Dirac-7 matrices are related to their 
Minkowski counterparts by 

7 4 = 74 = 7m, 7fe = l k = -ilM, (11) 
{ 7m ,7„} = 2V, 7^ = 7^. (12) 
7 5 = 75 = 74717273 = 7m- (13) 

Throughout this paper, we use the standard Dirac-Pauli 
representation for the 7-matrices: 





■° 








lk = ! 




l 


74 = 


0" 



where the Pauli spin matrices are given by 

= (? J)' ff2 = (i ~o)' ct3 = (0 -!)■ (l5) 

It has long been known that operators constructed out 
of smeared fields have dramatically reduced mixings with 
the high frequency modes of the theory. Thus, our oper- 
ators are constructed using spatially-smoothed link vari- 
ables Uj{x) and spatially-smeared quark fields ^(a;). The 
spatial links are smeared using either the stout-link pro- 
cedure described in Ref. [23 or the method introduced in 
Ref. j2^|. Note that only spatial staples are used in the 
link smoothening; no temporal staples are used, and the 
temporal link variables are not smeared. The smeared 
quark fields are defined bvpij 

$(x) = (l + A) ° V(*0, (16) 
1>(x) = m \l + -JL- A J , (17) 
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where a s and n„ are tunable parameters (n„ is a posi- 
tive integer) and the three-dimensional covariant Lapla- 
cian operators are defined in terms of the smeared link 
variables Uj(x) as follows: 

AO(s) = Yl (Uk(x)0(x+k) - 0(x)\ (18) 

ft=±l,±2,±3^ ' 

D»A = J] (0(^+^)^(^-0(^,(19) 

fc=±l,±2,±3^ ' 

where Oix) and 0(x) are operators defined at lattice 
site x with appropriate color structure, and noting that 

U-k{x) = U\(x — k). The smeared fields -0 and ip are 
Grassmann- valued; in particular, these fields anticom- 
mute in the same way that the original fields do, and the 
square of each smeared field vanishes. 

Our operators are designed with one eye on capturing 
the states of interest and the other eye on facilitating the 
efficient computation of the hadron correlation functions. 
Since the baryon resonances are expected to be large ob- 
jects, the use of extended operators is crucial. Our choice 
of basic building blocks is motivated by the need to incor- 
porate spatial extensions so as to facilitate the efficient 
and gauge-invariant assembly of many large hadron op- 
erators. To capture orbital structure, the quarks must 
be displaced in different directions, and to capture radial 
structure, the quarks must be displaced by different dis- 
tances. To maintain gauge invariance, covariant displace- 
ments in terms of the gauge-field parallel transporters 
must be used. To simplify matters, we consider only 
straight-path displacements in the six directions along 
the axes of the spatial cubic lattice: j = ±1,±2, ±3. 
Hence, we shall assemble our baryon (and later, meson) 
operators using the following basic building blocks: 

t a , vol, $Df*)L> ; = ±i.±2.±3, 

where A is a flavor index, a is a color index, a is a Dirac 
spin index, and the p-link gauge-covariant displacement 
operator in the j-th direction is defined by 

Df\x,7>) = Ujix) Uj(x+j) . . . Ujix+ip-mS^^, 

(20) 

for j = ±1,±2,±3 and p > 1. In what follows, we can 
achieve a significant economy in the notation by defining 
a zero-displacement operator 

D$'\x,x')=5 xx ,, (21) 

to indicate no displacement. In this way, our basic build- 
ing blocks can be listed more succinctly: 

(?3Wt£, -3<i<3. (22) 

Sometimes it will be convenient to make the flavor quan- 
tum number more apparent by writing 

u aa (x)=^ a (x), d aa (x)=ifi a (x), (23) 



and similarly for the s, c, b quarks. It is important to 
remember that the u, d, H, etc. operators so defined refer 
to smeared quark fields. 

IV. THREE-QUARK ELEMENTAL 
OPERATORS 

Having chosen the basic building blocks for our hadron 
operators listed in Eq. I|22|l . the next step is to devise 
elemental baryon operators Bf (t) having the appropri- 
ate color and flavor structure. These operators are cho- 
sen such that they are gauge-invariant and transform ir- 
rcducibly under the isotopic spin symmetry. Explicitly 
dealing with SU(3) flavor symmetry is not necessary, as 
will be discussed below. 

Gauge invariance is easily handled. Given three quarks 
with color indices a, 6, c associated with the same lattice 
site x, there is only one way of combining the color indices 
to arrive at a locally gauge-invariant object: the use of 
the antisymmetric Levi-Civita symbol e a bc- Similarly, 
covariantly-displaced quark fields must also be connected 
by an e a bc coupling at a single lattice site. To simplify the 
operator construction, we consider only combinations of 
displaced quarks having the same displacement length p. 
Thus, all of our three-quark baryon operators are linear 
superpositions of gauge-invariant terms of the form 

= e abc ((51 p) v^ q om) 

x&M$)*(x,t) (24) 
The "barred" operators have the form 

Note that the barred composite operator $ is defined dif- 
ferently than the barred fermion field operator ip due to 
the presence of the 74 spin matrices in Eq. I|25[) . Their 
presence is needed to ensure that the resulting correla- 
tion matrices satisfy the desirable property of hermiticity. 
These 74 matrices do not affect the transformation prop- 
erties of these operators under proper spatial rotations 
and parity, although they do affect the Lorentz boost 
properties. 

The simplest way to combine the basic building blocks 
previously described is to combine three quark fields at 
a single lattice site, corresponding to i = j = k = in 
Eq. lj2"4l . We refer to these operators as single-site opera- 
tors. Next, one of the three quarks can be displaced; such 
singly-displaced operators correspond to i = j = 0, k ^ 
in Eq. (|24p. In these operators, the two quarks which 
are not displaced may be viewed as forming a localized 
diquark, so such operators may be important if baryon 
formation is dominated by a quark-diquark mechanism. 
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TABLE I: The six types of three-quark & a pS ijk operators 
used, where A, B, C indicate the quark flavors, 1 < a, f3, j < 4 
are Dirac spin indices, and —3 < i,j,k < 3 are displace- 
ment indices. Smeared quark fields are shown by solid circles, 
line segments indicate covariant displacements, and each hol- 
low circle indicates the location of a color e a bc coupling. To 
simplify matters, all displacements have the same length (p 
lattice links) in any given operator. Remember that a dis- 
placement index having a zero value indicates no displace- 
ment. The singly-displaced operators are meant to mock up 
a quark-diquark combination, and the doubly-displaced and 
triply-displaced operators are chosen since they may favor A- 
flux and F-flux configurations, respectively. 



Operator type 



Displacement indices 



® 

single-site 

(D-» 

singly-displaced 
doubly-displaced-I 



doubly-displaced-L 



triply-displaced-T 



i = j — k = 



i = j = 0, k / 



i = 0, j = -k, k / 



i = 0, \j\ jt \k\, jk^O 



i = -3, \j\ + |fc|, Jk^0 



\i\ / \j\ ^ \k\, ijk^O 



triply-displaced-0 



Two of the quarks can be displaced; they can be dis- 
placed either in opposite directions (doubly-displaced-I), 
or in orthogonal directions (doubly-displaced-L). In such 
operators, one can in general choose different lengths for 
the two displacements, but to simplify matters, we re- 
strict our attention to the case in which both displace- 
ments have the same length p. Since the displacement of 
two quark fields essentially results in an object in which 
all three quarks are at different lattice sites, one may 
be tempted to exclude these operators in favor of triply- 
displaced operators. However, the relative importance 
of so-called F-flux and A-flux formation of the gluon 
field in three-quark systems (the A is actually a quan- 
tum superposition of F-flux forms) is a much-discussed 
issue (see, for example, Refs. [H Hi M, HI E El ) . 
Hence, it is important to include some operators with 
significant overlap with a A-flux configuration, as well 
as operators with strong mixing with a K-flux configura- 
tion. The doubly-displaced operators above are suitable 
for A-flux formation. 

Lastly, all three quarks can be displaced (again, re- 
stricting attention to the case of equal distances). If two 
of the quarks are displaced in opposite directions, this 
produces a coplanar T-shape (triply-displaced-T); alter- 



TABLE II: Elemental baryonic operators which annihilate 
certain states of definite isospin I, maximal ^3 = /, and 
strangeness S, in terms of the gauge-invariant three-quark 
operators defined in Eq. 1241 . 



Baryon 



I = h 



S 



Annihilation operators 



A++ 


3 
2 





^aj3 r y; ijk 


£+ 


1 


-1 


ijk 


N+ 


1 

2 





^aft'y; ijk ^a{3-y; ijk 


cO 


1 

2 


-2 


^ap'y; ijk 


A" 





-1 


jf.uds fodus 

ijk ijk 







-3 


<b sss 

ijk 



natively, the three quarks can be displaced in mutually 
orthogonal directions (triply-displaced-O). These opera- 
tors are suitable for producing K-flux configurations. 

These operators, summarized and illustrated in Ta- 
ble allow a large number of baryon operators to be 
constructed using a relatively small number of quark 
propagator sources. For a given reference source site x, 
quark propagators must be evaluated using only a hand- 
ful of different sources: for each quark mass value, we 
need a local source and displaced sources in at most each 
of the six directions. However, rotational invariance of 
the baryon correlation functions can be exploited to re- 
duce the number of source displacement directions. For 
singly-displaced operators at the source, a simultaneous 
rotation of the source and sink can always be used to 
align the source displacement along the +z direction, say. 
For doubly-displaced-I sources, a rotation can always be 
done to align one displacement along the +z direction 
and the other along the — z direction. Doubly-displaccd- 
L sources can be rotated so the source displacements 
are along the +y and +z directions, triply-displaced-T 
sources can be rotated to align the displacements along 
the +y, +z, and — z directions, and triply-displaced-0 
sources can always be rotated so the displacements are 
along the +x, +y, +z directions. In total, only source 
displacements along the +x, +y, +z, — z directions (four 
directions) are required. Hence, the number of conjugate 
gradient inversions needed is 



inversions 



= iV c N sp N K (l + 4N p ), 



(26) 



where N c — 3 is the number of colors, N sp — 4 is the 
number of Dirac spin components, N p is the number of 
displacement lengths p, and N K is the number of quark 
masses to be used. 

Incorporating the isospin symmetry is also straightfor- 
wardly done. Let 7i,T2,T3 denote the three hermitian 
generators of the isospin symmetry satisfying the com- 
mutation relations [7$, Tj] = iSiju^k- A set of quantum 

operators O^/J transforms under isospin according to the 
irreducible representation / if 



[t 3 , O l3 



- T 77 (/) 



(27) 
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TABLE III: The A ++ elemental operators contain three u 
quarks and are given by $^3". as defined in Eq. 12411 with 
spin components 1 < a, (3, 7 < 4 and displacement indices 
—3 < i, j, k < 3. The linearly independent operators we chose 
are described in the second column, and the numbers of inde- 
pendent operators of each type are listed in the third column. 
Replacing each u quark by an s quark in the operators below 
yields fl~ elemental operators. 



Operator type 



Restrictions 



Multiplicity 



single-site 

singly-displaced 

doubly-displaced-I 

doubly- displaced- L 

triply-displaced-T 

triply-displaced-0 



(i=j = k = 0, a < (3 < 7) 20 

(z=j = 0,fc^0, a< p) 240 

(i = 0, j = -k, k> 0) 192 

(i = 0,7fc^0; |j|<|fc|) 768 

(i = -j, j >0, 3 + |fc|, 768 

(*j*/0, |»|<|j|<|*|) 512 



[t+, Ol>] = v /(7-7 3 )(7 + 7 3 + l) (28) 



— (i) 
H 1 ) 



[t-, Ol>] = v/(J + 7 3 )(/-/3 + l) 0l ] _ x , (29) 
where t± = 17 ± 1T2 . It follows from these relations that 

[rs, [t 3 , 0^]] + i[r + , [r_, O^]] 

+i[r_, [r + , 0lf]] = 7(7+1) o£>. (30) 

Baryonic operators of definite isospin 7 and 7 3 are easily 
constructed using the above relations and the following 
commutation rules involving the isospin generators and 
the barred u, d, s quark field operators (suppressing all 
indices except flavor): 

[r 3 , u] = |u, [T3, d] = -fd, [t 3 , s] = 0, 
[t+,5] = 0, [r+, d] = u, [t+,S] = 0, (31) 
[r_, u] = d, [r_, 3] = 0, [r_, s] = 0, 



TABLE V: The N + elemental destruction operators are given 
(see Eq. I|24[l i. The linearly indepen- 



by $r d 



duu 
a/37; ijk 



dent operators we chose are described in the second column, 
and the numbers of independent operators of each type are 
listed in the third column. 



Operator type 



Restrictions 



Multiplicity 



single-site 

singly-displaced 

doubly-displaced-I 

doubly-displaced-L 

triply-displaced-T 



=j = k = 0, a > p, a > 7) 20 

(i = j = 0,fc/0) 384 

(i = 0, j — —k, k± 0) 384 

(i = 0,jk^0; |j| + |fc|) 1536 

[i = -j, jk^O, \j\^\k\) 1536 



triply-displaced-O (ijk^=0, \i\<\j\, |i|<|fc|, \j\ / |fc|) 1024 



and for the unbarred field operators, 

[t 3 , u] = -\u, [t 3 , d] = |d, [r 3 , s] = 0, 

[r_, u] = 0, [r_, d] = [r_, s] = 0, (32) 

[r + , u] = -d, [t+, d] = 0, [r+, s] = 0. 

Due to the isospin symmetry of QCD in the m u = m<j 
approximation, the particle masses do not depend on 7 3 , 
so we limit our attention to only one value of 7 3 in each 
isospin 7, strangeness S sector, choosing the maximal 
7 3 = 7 value. Using the above isospin relations, we first 
write down all possible flavor combinations appropriate 
for each isospin channel for the three-quark elemental op- 
erators from Table ID These are listed in Table [H] For 
each flavor sector and quark-displacement type, we then 
determine a maximal set of linearly independent opera- 
tors, giving us the final set of elemental operators we use. 
A symbolic computer program capable of manipulating 
Grassmann fields was written using Maple 9.5 and uti- 
lized to identify linearly independent operators. The in- 
dependent elemental operators we chose in the different 
flavor sections are described in Tables ITTT1 II VI Ivl and I VII 



TABLE IV: The E + elemental operators are given by 3 > ™^ ; yfc 
as defined in Eq. 12 H . The linearly independent operators we 
chose are described in the second column, and the numbers 
of independent operators of each type are listed in the third 
column. Interchanging the u and s quarks yields the H° ele- 
mental operators. 



Operator type 


Restrictions Multiplicity 


single-site 


(i=j = k = 0, a<P) 


40 


singly-displaced 


(i=j = 0,k=£0,a</3) 


240 


singly-displaced 


(i = k = 0, jjtO) 


384 


doubly-displaced-I 


(i = 0,j = -k,k?0) 


384 


doubly-displaced-I 


{i = -j, j > 0, *=0) 


192 


doubly-displaced-L 


(i = 0, jk^O, \j\ j4 |ft|) 


1536 


doubly-displaced-L 


(v /o, |i| < lil, fc=o) 


768 


triply-displaced-T 


{i=-j,j>0,j^\k\,k^0) 


768 


triply-displaced-T 


(ii/0, \i\^\j\,k = -i) 


1536 


triply-displaced-O 


(»ifc^o,|i|<|i|, K|^|fe|, \Mk\) 


1536 



TABLE VI: The A elemental destruction operators are given 
by ®a%;ijk - ®%g r ,ijk (see Eq. ED)- The linearly indepen- 
dent operators we chose are described in the second column, 
and the numbers of independent operators of each type are 
listed in the third column. 

Operator type Restrictions Multiplicity 



single-site (i=j = k = Q, a < 0) 24 

singly-displaced (i = j = 0, k 7^ 0, a < 0) 144 

singly-displaced (i — k — 0, j 7^ 0) 384 

doubly-displaced-I (i=0, j = -k, k ^ 0) 384 

doubly-displaced-I (i=-j, j > 0, k=0) 192 

doubly-displaced-L {i~0, jk 7^ 0, \j\ 7^ |*|) 1536 

doubly-displaced-L (ij 7^ 0, |i| < \ j\, fe = 0) 768 

triply-displaced-T (i = -j, j > 0, j ^ |fc|, k 0) 768 

triply-displaced-T (ij / 0, \i\ / \j\, k = -i) 1536 

triply-displaced-O (ijk^O, \i\<\j\, 1536 
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TABLE VII: Continuum limit spin identification: the number 
n J K of times that the A irrep of the octahedral point group 
Oh occurs in the (reducible) subduction of the J irrep of 
5(7(2). The numbers for G\ u ,G2u, H u are the same as for 
Gig , Gig , H g , respectively. 



J 




n G 2a 


J 


J 


n G lo 




n a 


1 

2 


1 








9 
2 


1 





2 


3 
2 








1 


11 

2 


1 


1 


2 


5 
2 





1 


1 


13 
2 


1 


2 


2 


7 
2 


1 


1 


1 


IS 
2 


1 


1 


3 



Due to an approximate SU (3) itds-flavor symmetry, 
quark flavor combinations in baryon operators are usu- 
ally chosen according to the irreducible representations 
of SU(3) flavor. Such combinations are simply linear 
superpositions of the operators presented above. Since 
we plan to obtain Monte Carlo estimates of the complete 
correlation matrices of operators including all allowed fla- 
vor combinations, the use of linear superpositions which 
transform irreducibly under SU (3) flavor is unnecessary. 
Our choice of operators described above is dictated by 
computational simplicity and efficiency. 

Lastly, note that the baryons A c , S c , S cc , and fl ccc can 
be studied using the operators presented above if the s 
quark is replaced with a c quark. Similarly, replacing the 
s quark by a b quark in the above operators allows us to 
investigate the Ah, E&, Sbb, and tlhbb baryons. Some other 
baryons, such as Q, cc , can be studied using other suit- 
able flavor replacements in the above operators, whereas 
the investigations of baryons such as S c containing use 
quarks cannot directly exploit the above tables, requiring 
slight modifications. 

V. PROJECTIONS ONTO SYMMETRY 
SECTORS 

The final step in our operator construction is to apply 
group-theoretical projections to obtain operators which 
transform irreducibly under all lattice rotation and re- 
flection symmetries. The basic building blocks used to 
assemble our baryon operators transform under the al- 
lowed spatial rotations and reflections of the point group 
Oh according to 

u R (D?^( X )) A y R = s{R) -i (scg^jto))^ (33) 

U R ^(x)D^Xa U R = ^)5j )t )i5(%,(34) 

where the transformation matrices for spatial inversion 
I s and proper rotations C n j through angle 2n/n about 
axis Oj are given by 

s ( c nj) = exp(hw tiV ['y ti ,'y v ]) ) (35) 
S(I S ) = 74, (36) 



with Uki = —2TtEjki/n and u>4k = t^ki = (lu^ u is an 
antisymmetric tensor which parametrizes rotations and 
boosts). A rotation by 7r/2 about the y-axis is conven- 
tionally denoted by C^y, and C^ z denotes a rotation by 
it/2 about the z-axis. These particular group elements 
are given by 

S(C iy ) = -L(l +7173 ), S(C iz ) = -L(l+ 727l ). (37) 

The allowed rotations on a three-dimensional spatially- 
isotropic cubic lattice form the octahedral group O which 
has 24 elements. Inclusion of spatial inversion yields the 
point group Oh which has 48 elements occurring in ten 
conjugacy classes. All elements of Oh can be generated 
from appropriate products of only C^y, Ci Zl and I s . 

Charge conjugation is another symmetry of our theory. 
Under charge conjugation C, the link variables U — > U* 
and our basic building blocks transform according to 

C(5f^))l C t = mDf%Cl, (38) 

where the charge conjugation matrix C must be unitary, 
antisymmetric, and satisfy C 7/J C^ = — 7 J. Our choice 
for C in the Dirac-Pauli representation is C — 7472- 

Operators which transform according to the irreducible 
representations of Oh can then be constructed using 
the well-known group-theoretical projections given in 
Eq. (jSJ. Orthogonality relations and hence, projection 
techniques, in group theory apply only to single- valued ir- 
reducible representations. However, the fermionic repre- 
sentations are double- valued representations of Oh- The 
commonly-used trick to circumvent this difficulty is to 
exploit the equivalence of the double- valued irreps of Oh 
with the extra single- valued irreps of the so-called double 
point group Off. This group is formed by introducing a 

new element E which represents a rotation by an angle 

2 

2tt about any axis, such that E = E (the identity). By 
including such an element, the total number of elements 
in Off is double the number of elements in Oh- The 96 
elements of Off occur in sixteen conjugacy classes. 

Since baryons are fermions, we need only be concerned 
with the six double- valued irreps of Oh- There are four 
two-dimensional irreps G\ g , G\ u , Gi g , and G2U, and two 
four-dimensional irreps H g and H u . The subscript g 
refers to even parity states, whereas the subscript u refers 
to odd parity states. The irreps G\ g and G± u contain the 
spin-1/2 states, spin-3/2 states reside in the H g and H u , 
and two of the spin projections of the spin-5/2 states oc- 
cur in the G 2ff and Giu irreps, while the remaining four 
projections reside in the H g and H u irreps. The spin 
content of each Oh irrep in the continuum limit is sum- 
marized in Table lvTTl This table lists the number of times 
that each of the Oh irreps occurs in the J = |, | , |, • • • 
representations of SU(2) subduced to Oh- 

To carry out the projections in Eq. ©, explicit rep- 
resentation matrices are needed. Our choice of represen- 
tation matrices is summarized in Table I Villi Matrices 
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TABLE VIII: Our choice of the representation matrices for the 
double- valued irreps of Oh- The Gi u , G211, H u matrices for the 
rotations Ciy, d z are the same as the Gi g , Gig, H g matrices, 
respectively, given below. Each of the Gi g ,G2 g , H g matrices 
for spatial inversion I s is the identity matrix, whereas each 
of the Giu,G2u, H u matrices for I 3 is — 1 times the identity 
matrix. The matrices for all other group elements can be ob- 
tained from appropriate multiplications of the Ci y , Ci z , and 
I s matrices. 

A r {A) (c iy ) r w {c iz ) 





1 


' 1 -1 " 






1 


"l-i 






71 


1 1 _ 








1+i _ 






-1 


' 1 -1 " 






-1 


' l-i 






71 


1 1 _ 






71 


1+i _ 






1 - 


V3 ^3 


-1 " 




"-1- 


-i 


1 


Vs 


-1 -1 


V3 


1 





l-i 






1 -1 


-y/3 


71 





1+i 




1 


^3 V$ 


1 . 











for only the group elements C±y, Ci Z , and I s are given in 
Table IVlTll since the representation matrices for all other 
group elements can be obtained by suitable multiplica- 
tions of the matrices for the three generating elements. 
For baryons, the representation matrix for E in each of 
the extra irreps is —1 times the identity matrix. 

Note that Table IVIII is the key to identifying the 
continuum-limit spin J corresponding to the masses ex- 
tracted in our Monte Carlo calculations. For example, to 
identify an even parity baryon as having J = 1/2, a level 
must be observed in the G\ g channel, and there must be 
no degenerate partners in either of the Gi g or H g . A level 
observed in the H g channel with no degenerate partners 
in the G\ g and G2 g channels (in the continuum limit) is 
a J = 3/2 state. Degenerate partners observed in the 
Gi g and H g channels with no partner in the G\ g channel 
indicates a J = 5/2 baryon. In other words, Table IVTll 
details the patterns of continuum-limit degeneracies cor- 
responding to each half-integral J value. 

Our operators are constructed using fcrmion fields 
ip(x) which annihilate a quark and create an antiquark. 
Hence, each of our baryon operators annihilates a three- 
quark system of a given parity P and creates a three- 
antiquark system of the same parity P. This means that 
in the baryon propagator, a baryon of parity P propa- 
gates forward in time while an antibaryon of parity P 
propagates backwards in time. Unlike boson fields, a 
fcrmion and its antifermion have opposite intrinsic par- 
ities, so that the antibaryon propagating backwards in 
time is the antiparticle of the parity partner of the baryon 
propagating forwards in time. Since chiral symmetry is 
spontaneously broken, the masses of parity partners may 
differ. The forward propagating baryon will have a mass 
different from that of the antibaryon propagating back- 
wards in time. If the even and odd parity baryon opera- 



tors are carefully designed with respect to one another, it 
is possible to arrange a definite relationship between the 
correlation matrix elements of one parity for t > and 
the opposite-parity matrix elements for t < 0, allowing 
an increase in statistics. Our operators are designed to 
take advantage of this symmetry (see below) . 

Our construction of the irreducible Bf XF (t) baryon op- 
erators is done in the following sequence of steps. 

(a) A set of Mb linearly-independent elemental opera- 
tors Bj (t) that transform among one another under 

is identified with the help of the computer software men- 
tioned earlier. This is done by starting with all possible 
operators of a given type (single-site, singly-displaced, 
and so on), then using the Maple program to detect 
dependencies between the operators. To find such de- 
pendencies, the computer program expresses each Bf (t) 
operator as a sum of products of Grassmann fields (or 
gauge-covariantly displaced Grassmann fields) with ex- 
plicit color, flavor, and spin indices: Bf = J2k dfk^k- 
The coefficients gf k form a Mb x Mo matrix, where Mo 
is the total number of Ok operators encountered, which is 
much larger than Mb ■ Since each Ok operator is a single 
product of three displaced Grassmann fields, these opera- 
tors are linearly independent, so the linear independence 
of the Bf operators boils down to the linear indepen- 
dence of the rows of the gf k matrix. In this way, the set 
of all possible operators of a given type is easily reduced 
to a set containing only linearly- independent operators. 

(b) We obtain the Mb x Mb representation matrices 
Wij (R) which satisfy 

Mb 

U R Bf{t)ul = Y J Bf(t)W Jl (R). (40) 
3=1 

Our Maple program determines the i-th column of the 
Wji(R) matrix for a given symmetry transformation R 
as follows. First, the Bf (t) operator is expressed as a 
sum of products of displaced Grassmann fields with ex- 
plicit color, flavor, and spin indices as in the previous 
step: Bf — J2kdfkOk- Next, the R symmetry trans- 
formation is applied to the displaced Grassmann fields 
in each Ok term in this sum of products using Eq. H33|) . 
The resulting sum of terms Ur Bf(t) U R — J^k^fkOk 
is then expressed as a linear superposition of the original 
Mb operators using the Moore-Penrose pseudoinverse [32^ 
of the gf k matrix. If the transformed operator contains 
Grassmann products which are not in the original set of 
Ok operators, then this signals that the starting basis of 
Bj operators is incomplete, but with our method in the 
previous step of choosing linearly-independent operators, 
this did not occur. In this way, the Wji(R) matrices for 
the generating group elements C± v , Ci z , and I s are ob- 
tained. The matrices for all other group elements are 
then determined by appropriate products of these three 
matrices. At this point, we have a set of Mb operators 
Bf (t) which form the basis of a reducible representation 
given by the Wij(R) matrices. Our remaining task is to 
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find a change of basis such that the resulting represen- 
tation matrices are block-diagonal with the blocks given 
by the irreducible representations of Oh- 

(c) Since the resulting Wji(R) matrices may not be 
unitary, we compute the hermitian metric matrix M 



Mi. 



^ M B 

= — E Z)Wfci(ii)*w fci (ii). 



So? 



(41) 



ReoP k=l 



This matrix will be needed in a later step where it will 
facilitate the full block-diagonalization of the Wij (R) ma- 
trices. 

(d) For each even-parity irrep A, we compute the large 
Mb x Mb projection matrix for row A = 1: 



pKXF _ 



9oP 



Reot 



[T^(R)W ji (R) 



A=l 



This is one of the most important steps in our operator 
construction. Applying the group-theoretical projection 
of Eq. © to the operator Bf, then utilizing Eq. |@DJ|, 
produces a new operator Bf XF {t) = J^j Pij Xf (*) 
which resides in the subspace of operators which trans- 
form according to the row A = 1 of the given irrep A. 
In other words, the i-th row of the projection matrix P 
contains the superposition coefficients of the projected 
gAAF p era tor. 

(e) Although group theory guarantees that the result- 
ing projected Bf XF operators reside in the subspace as- 
sociated with row A of the A irrep, it does not guarantee 
that all of the resulting operators are linearly indepen- 
dent. The maximum number of independent operators 
in the projected subspace is given by the rank r of the 
projection matrix. Hence, the next step is to form r 
superpositions of the Bf XF operators such that the re- 
sulting r operators are linearly independent. The choice 
of these operators is not unique. In practice, these linear 
combinations are obtained using the well-known Gram- 
Schmidt procedure, but with a modified inner product to 
incorporate the metric matrix M. The use of the metric 
matrix M ensures full block-diagonalization of the orig- 
inal Wij(R) matrices. Using such a procedure, the final 
operators, expressed in terms of the original set of Bf 
operators by 



M B 



B^ F (t)=J2 & F B F (t), (A = l) 



(43) 



have superposition coefficients c^ XF that satisfy 

M B 



k,l=l 



AAF* 



Ma c 



kl c it 



AAF 



Sij, (i = l...r). 



(44) 



(f) For each of the r operators Bf XF (t) in the first 
row A = 1, we obtain partner operators in all other rows 



jti > 1 using 



M B 

£' 

3=1 



,AAF 



9oP 



Reot 



T^(R)W kj (R). 



(45) 



The use of operators belonging to other rows will be im- 
portant for increasing the statistics of our Monte Carlo 
calculations, as will be discussed below. 

(g) Although the odd-parity operators can be obtained 
using the same procedure described above, we instead 
utilize charge conjugation to construct the odd-parity op- 
erators. Consider the correlation matrix element of even- 
parity operators for t > 0. Suppressing flavor and dis- 
placement indices, one sees that invariance under charge 
conjugation implies that 



(42) Cy( t ) 



c« 7 cg:<0| $ Q , 7 (0$^(0) |0), 

c « 7 cg;<o| etc^ 7 (i)ctc%_(o)ctc | 0) 



= cg 7 cgl<0| ^(^(0) |0) 
= c« 7 cgl<0| $ a ,^ y (0)$^_,(-i) |0) 

2 2 2 2 2 2 
X ^a' a"/ /3' /3"/j' j^ta'a^f'p'"p"f : y'' : y^ 

using invariance under time translations of the above 
expectation value and invariance of the vacuum under 
charge conjugation. The last line above represents the 
correlation of odd-parity operators propagating tempo- 
rally backwards. Hence, for a given even-parity operator 
Bf(t), we can define an odd-parity operator Bf(t) by 
rotating the three Dirac indices using the 72 matrix and 
replacing the expansion coefficients by their complex con- 
jugates such that the correlation matrices of the even and 
odd parity operators are related by 



r 



(*) = -cJH-*)* 



(46) 



For a lattice of N t sites in the time direction with peri- 
odic {rj t = 1) or antiperiodic (rj t = —1) boundary condi- 
tions, this means that 



{N t - ty 



(47) 



and similarly for the other irreps. This allows us to ap- 
propriately average over forward and backward temporal 
propagations for increased statistics. 

In the absence of any external applied fields, the ener- 
gies of the baryons do not depend on the row A of a given 
irrep A, so we can increase statistics by averaging over all 
rows. The construction of the operators in the different 
rows of the irreps as described above leads to correlation 
matrices which satisfy Cfj XF (t) = Cf^ lF (t), for all rows 
A, fi. Hence, the correlation matrix elements themselves 
can be averaged over rows. 
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TABLE IX: The single-site A ++ operators which transform 
irreducibly under the symmetry group of the spatial lattice, 
defining A af3j = $^. 0Q0 (see Eq. ^). 



Irrep 


Row 


Operator 


Row 


Operator 


Gig 


1 


A134 — A233 


2 


A144 — A234 


Giu 


1 


A123 — A114 


2 


A223 — A124 


H g 


1 


A ln 


2 


\/3An2 


H g 


3 


\/3Ai 22 


4 


A222 


H g 


1 


Vs A133 


2 


2A134 + A233 


H a 


3 


A144 + 2A234 


4 


^3 A244 


H u 


1 


A333 


2 


n/3 A334 


Hu 


3 


\/3 A344 


4 


A444 


H u 


1 


V3An3 


2 


A114 + 2A123 


H v 


3 


2A124 + A223 


4 


n/3 A224 



The single-site operators produced by the above pro- 
cedure are presented in Tables IIXIXIII It is not possi- 
ble to list all of the singly-displaced, doubly-displaced, 
and triply-displaced operators in this paper. Instead, we 
simply list the numbers of operators of each type which 
project into each row of the irreps for the A ++ , S + , iV + , 
and A baryons in Table IXIIII Further details about 
these operators are available upon request. 

Note that we have not yet attempted to remove re- 
dundant operators related to others by a total lattice 
derivative, and that there may exist relationships be- 
tween the correlation matrix elements of these operators 
from the equations of motion. Such relationships can 
be easily identified using a singular value decomposition 
of the correlation matrices in small-lattice low-statistics 
Monte Carlo calculations. Since the Schwinger-Dyson 
equations may relate operators of different types, such 
as singly-displaced operators with single-site operators, 



TABLE X: The single-site A operators which transform ir- 
reducibly under the symmetry group of the spatial lattice, 
defining A Q/37 = ^%. 000 - $^000 (see Eq. $ZQ). 



Irrep 


Row 


Operator 


Row 


Operator 


Gi 9 


1 


A121 


2 


A122 


Gig 


1 


A341 


2 


A342 


Gig 


1 


A134 — A143 


2 


A234 — A243 


Gig 


1 


A134 +A143 — 2A233 


2 


2A144 — A234 — A243 


Giu 


1 


A343 


2 


A344 


Giu 


1 


A123 


2 


A124 


Giu 


1 


A231 — A132 


2 


A241 — A142 


Giu 


1 


2A141 — A132 — A231 


2 


A142 — 2A232 + A241 


H B 


1 


v^Aias 


2 


A134 + A143 + A233 


H g 


3 


A144+A2344-A243 


4 


y/3 A244 


H„ 


1 


-V3Aisi 


2 


— A132 — A141 — A23I 


H u 


3 


— A142 — A232 — A241 


4 


— \/3 A242 



TABLE XI: The single-site E + operators which transform ir- 
reducibly under the symmetry group of the spatial lattice, 
defining E a(37 = $^* ;00 o (see Eq. JStj). 



Irrep Row Operator Row Operator 

Gi s 

Gig 

Gig 
Gig 

Giu 
Giu 
Giu 
Giu 
H g 
H g 
H g 
H g 
H g 
H g 
H u 
H u 
H u 
H u 



identifying these relationships must be done at a later 
stage in the calculations. 



VI. BARYON PROPAGATORS 

To extract the baryon masses, we need to compute the 
correlations C± XF (t) = (0| TB^ XF {t)B j (0) |0), using 



TABLE XII: The single-site N + operators which transform 
irreducibly under the symmetry group of the spatial lattice, 



defining N a ^ = $Xooo - $Xooo (see Eq. 


Irrep 


Row 


Operator 


Row 


Operator 


Gig 


1 


N211 


2 


A221 


Gig 


1 


Nii 3 


2 


A 423 


Gig 


1 


2A 3 32 +A 4 i 3 -2Af 4 3i 


2 


2A 4 32-2A 44 l-A r 423 


Giu 


1 


A f 433 


2 


A443 


Giu 


1 


N321 — A312 


2 


A421 — N412 


Giu 


1 


A312 4-^321 -2AT 4U 


2 


2N 3 22 -N 4 i2-N 42 l 


Hg 


1 


V3N 3S i 


2 


A r 332-A 4 l 3 4-2A''431 


Hg 


3 


2^432 4^441-^423 


4 


%/3A 4 42 


H u 


1 


-v^Aan 


2 


-A 3 i2-7V 3 21-A i '411 


Hu 


3 


—N322 — N412 — N421 


4 


-^^422 



1 


5^112 — ^121 


2 


Sl22 — ^221 


1 


S134 — S143 


2 


Xi234 — ^243 


1 


S332 — S341 


2 


S342 — S441 


1 


Sl34-{-Sl43 — 2Xi233 


2 


2Sl44 — X/234 — 5j243 


1 


S334 — S343 


2 


S344 — S443 


1 


/ j1 91 y jOQI 

^lo/ ^zol 


2 




1 


S114 — £123 


2 


S124 — S223 


1 


S132 — 2E141 4-E231 


2 


2S232 — £241 — S142 


1 


\/3 Em 


2 


S112 4- 2E121 


3 


2E122 4- £221 


4 


V3 E222 


1 


v3 E133 


2 


Ei 3 4 4- S143 4-E233 


3 


E144 4- S234 4- S243 


4 


v3 E244 


1 


V% S331 


2 


E332 4- 2E341 


3 


2E342 4- E441 


4 


V 7 ^ E442 


1 


a/3 £333 


2 


E334 4- 2E343 


3 


2E344 4- S443 


4 


E444 


1 


V% £131 


2 


E132 4-E141 4-E231 


3 


E142 4-S232 4-S241 


4 


v3 E242 


1 


\/3 E113 


2 


E114 4- 2E123 


3 


2S124 4- E223 


4 


v 7 ^ E224 
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TABLE XIII: The numbers of operators of each type which project into each row of the Gi g , G?2 9 , and H g irreps for the A ++ , 
£ + , N + , and A baryons. The numbers for the Gi u , G2u, and H u irreps are the same as for the Gi g , G2 g , and H g , respectively. 







A++ 






£+ 






iV+ 






A 




Operator type 


Gi g 


G 2g 


H g 


Gi g 


G 2g 


Hg 


Gig 


G 2g 


H g 


Gig 


G 2 g 


H g 


single-site 


1 





2 


4 





3 


3 





1 


4 





1 


singly-displaced 


14 


6 


20 


38 


14 


52 


24 


8 


32 


34 


10 


44 


doubly-displaced-I 


12 


4 


16 


36 


12 


48 


24 


8 


32 


36 


12 


48 


doubly-displaced-L 


32 


32 


64 


96 


96 


192 


64 


64 


128 


96 


96 


192 


triply-displaced-T 


32 


32 


64 


96 


96 


192 


64 


64 


128 


96 


96 


192 


triply-displaced-0 


20 


20 


44 


64 


64 


128 


44 


44 


84 


64 


64 


128 



the operators constructed as described above. In this 
section, several issues related to computing these cor- 
relation matrix elements are discussed. In particular, 
we discuss the use of symmetry to minimize the number 
of quark-propagator sources, the use of gauge-invariant 
three-quark propagators as an intermediate step in the 
baryon propagator determinations, and detail the appli- 
cation of Wick's theorem. 

The baryon propagators may be expressed in terms of 
the correlations of the elemental operators by 



M B 

£■ 

k, 1=1 



AXFAXF* 



<0| Tflf (i)Bf (0) |0). (48) 



Note that correlations between operators in different rows 
of the same irrep vanish. Since the number of elemental 
operators is large and the quark propagators are rather 
expensive to compute, it is very important to use sym- 
metry to reduce the number of quark-propagator sources. 
Given the invariance of the vacuum and the unitarity of 
the symmetry transformation operators, we know that 



<0| TB^(t)B l (0) |0> 

= (0| T U R B[(t)U R U R B?(0)U R |0>, 



(49) 



= £ W klk (R)W V i(Rr(0\TB^(t)B v (0)\0) 



k',i'=i 



for any group element R of Oh- Hence, for each source 

F 

B l (0), we can choose a group element Ri such that we 
minimize the total number of source elemental operators 
which must be considered. For example, consider the 
singly-displaced operators. We can choose an Ri such 
that the displaced quark in the source is always displaced 
in the +z direction. Similarly, a group element Ri can 
always be chosen to rotate each of the other types of 
operators into a specific orientation. 

The coefficients c^ AF in the baryon operators involve 
only the Dirac spin components and the quark displace- 
ment directions and are independent of the color indices 
and spatial sites. Thus, in calculating the baryon corre- 
lators, it is convenient to first calculate gauge-invariant 
three-quark propagators in which all summations over 



color indices and spatial sites have been done. A three- 
quark propagator is defined by 

^{ABC)(pp) _ , , 

= y^ £ abc£-abc Q {A) - ,_T_(a;,£|a;o,0) 

aoc t aaip\aaip y ' ' 

x 

x C P iw^'^°' 0) ^Xi-ip (a; ' i|a;o ' 0) ' (50) 

where Q^ A \ . -_(x,t\x o ,0) denotes the propagator for a 

aatp\aaip K 1 ' 

single smeared quark field of flavor A from source site Xq 
at time t = to sink site x at time t. At the sink, a 
denotes color, a is the Dirac spin index, i is the displace- 
ment direction, and p is the displacement length, and 
similarly at the source for a, a, i, and p, respectively. No- 
tice that the three-quark propagator is symmetric under 
interchange of all indices associated with the same flavor. 
As usual, translation invariance is invoked at the source 
so that summation over spatial sites is done only at the 
sink. These three-quark propagators are computed for 
all possible values of the six Dirac spin indices. 

Each baryon correlator is simply a linear superposi- 
tion of elements of the three-quark propagators. These 
superposition coefficients are calculated as follows: first, 
the baryon operators at the source and sink are expressed 
in terms of the elemental operators; next, Wick's theo- 
rem is applied to express the correlator as a large sum 
of three-quark propagator components; finally, symme- 
try operations are applied to minimize the number of 
source orientations, and the results are averaged over the 
rows of the representations. C++ code was written to 
perform these computations, and the resulting superpo- 
sition coefficients are stored in computer files which are 
subsequently used as input to the Monte Carlo runs. 

Wick's theorem is an important part of expressing the 
baryon correlators in terms of the three-quark propaga- 
tors. To simplify the notation in the following, let the 
indices fj,,u, r each represent a Dirac spin index and a 
displacement direction, and suppress the displacement 
lengths. Define c$ r = 0^*^,7^,7^,7^,, then the ele- 
ments of the baryon correlation matrix in the A ++ chan- 
nel are given in terms of three-quark propagator compo- 
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nents (before source-minimizing rotations) by 
C {A) (t) - r W r (3) I n {uuu) (t) 

<^ij [ l > - C ^T C Ju^\^ r (r\Jl)^\V)( IJ ,\T)\ l ) 



G 



(uuu) 

(r|7T)(i/|r)(^|F) 



(uuu) 

(t|77)(i/|/I)(/x|t) 



«(*) 



„(«««) / , r (»nn) , 



+ G (rf^M7J)(M|57)W}- 



(51) 



The N + correlators are expressed in terms of components 
of three-quark propagators by 



C {N) (t) = c w c^L/g 



(uud) 

(^\JL)(v\V)(t\t) 



+ ^(p\V){u 


IP)(r|r) 




l^)(r|M) 




It) (r |7I) 


p^(uud) 

- <->r(t 




p^(uud) 








p^(uud) 










}■ 





(52) 



and for the E + and A channels, one finds 



Clf\t) = » 



( ^(uus) /.\ 



c 



(A) 



+ ^(™F)V|M)(T|r)(*)}' 



(53) 



s-y(uds) ^(uds) 



+ G 



(lids) 



r|r)}- 



(54) 



VII. CONCLUSION 



We plan to undertake a comprehensive study of the 
spectrum of QCD using Monte Carlo computations. Our 
first goal in this study is to calculate the masses of as 
many low- lying hadron resonances as possible. Success- 
fully extracting these masses will depend crucially on us- 
ing carefully designed spatially-extended hadronic oper- 
ators. In this first paper, the construction of three-quark 
A, AT, S, A, f2, 5 baryon operators using group-theoretical 
projections was detailed. The operators were assem- 
bled out of gauge-covariantly displaced quark fields and 



transform according to the irreducible representations 
of the symmetry group of a spatial simple cubic lat- 
tice. Single-site, singly-displaced, doubly-displaced, and 
triply-displaced three-quark operators were considered. 
The guiding principles in devising our operators were 
maximizing overlaps with the states of interest while min- 
imizing the number of quark-propagator sources. Identi- 
fying the continuum-limit spins J of the states was ad- 
dressed, and various issues related to computing the cor- 
relation matrix elements of the baryon operators were 
discussed. 

Due to the complexity of these calculations and the im- 
portance of providing checks on our final results, we have 
been pursuing two different approaches to constructing 
the baryon operators. An alternative method of build- 
ing the baryon operators based on Clebsch-Gordan tech- 
niques will be presented elsewhere [lq|. The construction 
of meson and multi-hadron operators will be described 
in subsequent papers. 

The Monte Carlo software to evaluate the correlation 
matrix elements of these baryons operators has been writ- 
ten and thoroughly tested using a large variety of checks, 
including comparison with known results in the case of a 
uniform constant background gauge field. This software 
uses the Chroma Software System for Lattice QCD [3^ 
with QDP++ and QMP, developed under the Scientific 
Discovery through Advanced Computing initiative of the 
U.S. Department of Energy. Preliminary results have al- 
ready been presented in Refs. Although a very 
large number of baryon operators have been devised, it 
is not our intent to evaluate correlation matrices using 
all of these operators. Such calculations would not be 
feasible. The next important step in our study is to 
remove dynamically-redundant and ineffective operators 
using low-statistics Monte Carlo calculations, with the 
goal of finding some reasonably small subsets of operators 
adequate for extracting the low- lying masses of interest. 
Such calculations are currently in progress. 
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